Cooling a mechanical resonator via coupling to a tunable double quantum dot 
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We study the cooling of a mechanical resonator (MR) that is capacitively coupled to a double 
quantum dot (DQD). The MR is cooled by the dynamical backaction induced by the capacitive 
coupling between the DQD and the MR. The DQD is excited by a microwave field and afterwards 
a tunneling event results in the decay of the excited state of the DQD. An important advantage 
of this system is that both the energy level splitting and the decay rate of the DQD can be well 
tuned by varying the gate voltage. We find that the steady average occupancy, below unity, of the 
MR can be achieved by changing both the decay rate of the excited state and the detuning between 
the transition frequency of the DQD and the microwave frequency, in analogy to the laser sideband 
cooling of an atom or trapped ion in atomic physics. Our results show that the cooling of the MR 
to the ground state is experimentally implementable. 
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I. INTRODUCTION 

Mechanical resonators (MRs) are currently attracting 
considerable interest because of their potential appli- 
cations to high-precision displacement detection? 1 - mass 
detection^ and quantum measurement, 3 Recent techni- 
cal advances allow the fabrication of a MR with both a 
high quality factor (Q-factor) and a sufficient high fre- 
quency, approaching 1 GHz 4^ Such a MR provides a 
good platform for exploring various quantum phenomena 
and for observing the quantum-to-classical transition in 
macroscopic objects^ Moreover, quantized MRs could 
be useful in quantum information science. Indeed, the 
quantized motion of buckling nanoscale bars has been 
proposed for implementing qubits^ i 10 ' 11 The generation 
of entangled states ; 12 ' 13 squeezed states^ and quantum 
nondemolition measurements 15 using MRs have also been 
studied. However, to prepare an ideal ground state, the 
basic requirement is to be able to cool the MR to a state 
with a mean phonon number (n) <C 1. 

Numerous experiments cool ing a single MR have been 
reported recently (see, e.g., [TMlSIlMl^IlOIsilH^) ■ In 
these experiments, the MR and a fixed micromirror 
form an optical cavity and the MR is cooled by ei- 
ther radiation-pressure-induced backaction or bolometric 
backaction. Experimental results show that a single MR 
can be cooled down from room temperature to an effec- 
tive temperature of the order of 0.1 K (Refs.[l0and[l7l) or 
10 K (Refs-HIii^). However, for a MR with frequency 
~ 20 MHz, a temperature lower than 1 mK is required 
in order to drive the MR to the quantum regime. Thus, 
more effective cooling methods are needed, in addition to 
increasing the oscillation frequency of the MR. Moreover, 
besides the classical and semiclassical analyses of cooling 
a single MR via dynamical backaction ) 20 ' 23 ' 24 some quan- 
tum theories have also been developed) 25 ' 26 ' 27 ' 28 For ex- 
ample, in Refs. 25 and 26, it is predicted that a MR can 
be cooled down to its ground state when the frequency 
ui m of the MR is either comparable to or larger than the 



optical cavity's resonance linewidth T. In laser cooling, 
this corresponds to the sideband cooling of a bound atom 
or a trapped ion, where the lowest occupancy attainable 
is given by (n) w r 2 /16w^ <C 1, indicating that the MR 
can be most of the time in its ground state. Recently, 
the resolved-sideband cooling of a MR has been realized 
by coupling the MR to an optical resonant system^ 

Besides optomechanical cooling, an alternative way 
would be to cool the MR by coupling it, via an electronic 
coupling, to an electronic system. This provides the ad- 
vantage of fabricating and integrating the electronic de- 
vice into a cryogenic system. In principle, the electronic 
cooling of a MR can be achieved by several means, in- 
cluding coupling the MR to: (1) an optical quantum 
dotj2°i (2) a Josephson-juction superconducting quantum 
devic o 31 ' 32 ' 33 ' 34 ' 35 (which behaves like a superconduct- 
ing artificial atom 3 ^), and (3) a one-dimensional trans- 
mission line^ 7 - Moreover, the experimental cooling of a 
MR, via coupling it to a superconductor single-electron 
transistor 38 or to an LC circuit, 39 has also been reported. 



A. Cooling a mechanical resonator coupled to a 
double quantum dot 

In this work we propose an approach to cool a MR 
by coupling it to an electronic system: a double quan- 
tum dot (DQD). Indeed, the whole system consists of a 
DQD and a MR. The MR, together with another static 
plate, forms a gate capacitor adjacent to the left dot (see 
Fig. [T|) . The oscillation of the MR will modulate the ef- 
fective capacitance of this capacitor. In this way, the MR 
can be strongly coupled to the DQD. 

The cooling mechanism can be understood as follows. 
Two localized states in the DQD (|1) and |2)), with en- 
ergy level splitting fitoo , are driven by a microwave field of 
frequency ix>d (Fig. [J) . Similar to the resolved-sideband 
cooling of a trapped ion, the DQD resonant transition 
frequency is modulated by the oscillation of the MR, and 



the absorption spectrum consists of a series of sidebands 
at frequencies (u>q — ju) m ), where j = ±1, ±2, .... When 
the energy level splitting Tlloq of the DQD is tuned to sat- 
isfy the lowest sideband condition, i.e., u) m = WQ—uid, the 
excitation of the DQD from the ground state |1) to the 
excited state |2) absorbs a photon of energy h(u>o — u> m ). 
The subsequent decay of the DQD via electron tunnel- 
ing emits a photon of energy Kloq. Hence, each scat- 
tering process carries away the MR's vibrational energy 
by hu> m , or reduces the MR's quantum number n by 1. 
This cooling process is described by the state transition 
|l)|n> -► |2)|n-l) -> |l)|n- 1) (see Fig. ^. A series of 
cycles of this process leads to the cooling of the MR. On 
the other hand, the reverse process 1 1) \n) — ► |2) \n + 1) — » 
|l)|n + 1), increasing the phonon number, is suppressed 
because it is off-resonance. Like in Sisyphus cooling, 40 a 
cycle in one direction induces cooling, while the cycle in 
the other direction produces heating (see Fig. [21). 



A similar approach of cooling a semiconductor beam 
by coupling it to an optical QD was proposed in Ref. [30. 
Comparing it to this study, our approach has the follow- 
ing potential advantages: (i) the cooling system can be 
fabricated more easily. The MR here is used as a part 
of the capacitor, which is easier to fabricate compared 
to the proposal in Ref. [30J that embeds an optical QD 
in a nanoscale beam; (ii) the decay rate T of the upper- 
state of the resonant system (DQD) is just the rate of the 
electron tunneling to the electrode, which is tunable by 
varying the gate voltage; (iii) with a DQD, it is easy to 
achieve the lowest sideband condition, i.e., u> m = u>o — Ud, 
for cooling the MR by changing the energy level splitting 
Uojq via the gate voltage. 

In typical transport experiment s 41 ' 42 with QDs, the 
tunneling rate T ranges from 10 kHz to 10 GHz, while 
the fundamental frequency ui m of the MR is of the order— 
of 100 MHz. The resolved-sideband cooling regime, i.e., 
Lu m 3> r, can be reached by tuning the tunneling rate 
r. In this regime, our results show that the steady aver- 
age phonon occupancy of the MR can be of the order of 
(T/uj m ) 2 <C 1, indicating that the MR can be cooled to 
the ground state. 

This paper is organized as follows. In Sec. II, we intro- 
duce a model of the coupled MR-DQD system and derive 
its effective Hamiltonian. In Sec. Ill, we derive the mas- 
ter equation of the coupled MR-DQD system and then 
eliminate the DQD's degrees of freedom to obtain the 
master equation for the reduced density matrix of the 
MR. With this master equation for the MR, we study, 
in Sec. IV, the cooling of the MR by considering the 
steady average occupancy of the MR. Moreover, we an- 
alyze the steady phonon occupancy of the MR in the 
resolved-sideband cooling regime. Section V summarizes 
our conclusions. Furthermore, in Appendices A and B, 
we show the derivations of the master equations. 
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FIG. 1: (Color online) (a) Schematic diagram of a DQD con- 
nected to an electron source S and drain D via tunneling bar- 
riers. An oscillating plate (the MR) and another static plate 
form a capacitor on the left dot, which provides a capacitive 
coupling between the oscillating MR and the DQD. The en- 
ergy level of each dot is tunable by varying the gate voltage 
V g applied to the dot through the capacitor, (b) Transport 
process of an electron through a DQD: First, an electron tun- 
nels from the source to the left dot, and then a microwave 
field drives it to the right dot. Finally, it tunnels to the drain 
on the right side. 



II. MECHANICAL RESONATOR COUPLED TO 
A DOUBLE QUANTUM DOT 

A. Model 

The circuit diagram of a MR coupled to a DQD is 
shown in Fig. [lja). The DQD is connected to two elec- 
trodes by tunneling barriers. The bias voltage across the 
DQD is set such that the chemical potential of the left 
electrode /j,l is higher than that of the right electrode 
/i/j and thus electrons can tunnel from the left electrode 
to the right one through the DQD. We assume that the 
DQD is in the Coulomb regime, such that at most a sin- 
gle electron is allowed in the DQD. The corresponding 
electron states of the DQD are, respectively, the vacuum 
state |0), one electron in the left dot |1), and one electron 
in the right dot |2). Here we consider the case where the 
hopping strength between the two dots is much smaller 
than the energy level splitting of the two dots. To excite 
the electron from the left dot to the right one, we apply 
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FIG. 2: (Color online) Schematic diagram of the cooling 
(heating) process in the coupled MR-DQD system. When 
the DQD is excited by a red-detuned microwave field, i.e., 
A = u>d— u)q < 0, the anti-Stokes process (|l)|n) — ► |2)|n — 1)) 
is resonantly enhanced. A subsequent decay from the ex- 
cited state |2) to the ground state |1) reduces the energy 
of the MR by one quanta. This cools the MR because the 
emitted (blue) photon has more energy than the absorbed 
(green) photon. Due to the off-resonance, the Stokes process 
(|l)|n) — » |2)|n + 1)) is suppressed. However, in the reso- 
nant case, the Stokes process dominates and the cycle heats 
the system. The blue (red) vertical downward arrow shows 
the cooling (heating) process, which decreases (increases) the 
phonon occupancy n of the MR. 



a microwave field to the DQD [see Fig. [T^b)] ■ 

As shown in Fig. [Ha), the capacitor adjacent to the 
left dot is formed by a static plate and a single MR with 
a gate voltage V g \ applied to it. Thus, the displacement 
of the MR from its equilibrium position will modulate 
the capacitance c g \{x). For typical experimental param- 
eters, the displacement x of the MR is much less than 
the equilibrium distance d between the two plates, i.e., 
x <C d. Hence, the capacitance can be approximately 
given by 



Cgl(x) ~ C 9l( 1 -^) 



(1) 



where C g \ is the capacitance at X = 0. For the har- 
monically oscillating MR discussed here, the quantized 
displacement operator x can be written as (here % = 1) 



V Zmuj m 



(2) 



where m is the effective mass of the MR and b + (b) is 
the bosonic creation (annihilation) operator. The total 
Hamiltonian of the whole system reads 



H, 



total : 



H + Hi n t + Ht, 



where 



H = -fflcads + #DQD + Hr, 



(3) 



(4) 



is the sum of the isolated bath Hamiltonian i?i oa ds, the 
Hamiltonian i^DQD of the DQD driven by a microwave 
field, and the Hamiltonian Hr of the MR, with 



-^ leads / -£ / qfcC" a fcC"qfc; 



Hbqd 



Wo 



a z + Qa x + Slo cos(wdi)iXz 



(5) 

(6) 

(7) 



Here 



J 



(Caka) is the creation (annihilation) operator 



of an electron with momentum k in electrode a (a — I, r). 



a\a 2 



a\ai and a x 



a 2 ai 



a[a2 are the Pauli 



matrices with a\ (a 2 ) being the electron creation operator 
in the left (right) dot of the DQD. The second term in 
Eq. ([5]), Q<J X i is the hopping tunneling term between the 
two dots. The third term in Eq. ([6]) describes the applied 
microwave field with driving frequency cud and amplitude 
Q . 

The coupling between the MR and the electron in the 
left dot is given by 45 



(8) 



H ln t = -A a\ ai (fe T + 6), 

with an electromechanical coupling strength A = r/uj m . 
For a typical electromechanical coupling, r\ ~ 10 _1 . The 
tunneling coupling between the DQD and the electrodes 

is 



H T = ^{fllk a\ Qfc + flrk a\ Crk + H.c), 



(9) 



where r2; fc ( rfc \ characterizes the coupling strength be- 
tween the QD and the left (right) lead. Hereafter, the 
subscript "Z" ("r") refers to the left (right) electrode. 



B. Effective Hamiltonian 

It is difficult to directly analyze the coupled system 
due to the different time scales for the dynamics of the 
DQD, the MR, and the coupling between them. In or- 
der to solve this problem, we first eliminate the coupling 
term between the MR and the DQD through a canonical 
transform U — e s on the whole system, where 

S = exp[— ?]a[ai(b' — b)]. 
With the following relations: 

UaxU* = ai exp[ri(tf -b)], 

Ua 2 U^ = a 2 , 
UbU^ = b + r)a\ai, 
the transformed Hamiltonian is given by 
H = UHtot^ 



(10) 



(11) 



E ak c' ak c ak + —a 



M 



-[Q + tt cos{uj d t)](cr+B^ + H.c.) 

- y^(Q; fe a^cikB^ + Q rk a\ c rk + H.c), (12) 



where we have redefined the energy level splitting ojq and 
5 = cxp[r/O t -&)]. (13) 



Here we introduce the ladder operators a + = a 2 a\ and 
<7_ = crjj_. 

Moreover, to eliminate the time-dependence of the 
driving term in Eq. (112|) , we now employ a unitary trans- 
form 



Ur = exp 






(14) 



to change the Hamiltonian to a rotating frame. By ne- 
glecting the fast-oscillating terms within the rotating- 
wave approximation, the resulting Hamiltonian is given 
by 

H = H sys + 2_^ E a kC ak c a k 

a k 

+ } j [Clik a+cikB^ exp (i-^-t) + ^rfe a\crk + H.c], 



(15) 



with 



#s 



rfb - 



^(a + St+H.c). (16) 



Here, A = uj d — ljq is the driving frequency detuning 
the microwave-field from the transition frequency of the 
DQD. The rotating-wave approximation is valid when 
LUd 3> f2, which corresponds to a weak hopping tunneling 
between the two dots. 



III. MASTER EQUATION FOR THE 
MECHANICAL RESONATOR 

Within the Born-Markov approximation, by integrat- 
ing over the electrode degrees of freedom, we derive a 
master equation for the coupled MR-DQD system 



dp 
~db 



-i[H sys , p] + £ T p + Cup- 



Here the Liouvillian operator Ct presents the tunneling 
events through the DQD in the presence of a single MR. 
In Eq. dnj) and below, the subscript "T" ("D") refers to 
the tunneling (dissipation). By expanding B in Eq. (JT3J) 
up to second order in 77 and assuming that the energy lev- 
els with one-phonon-mediated tunneling arc within the 
bias window^ CtP reads 

C TP = T^l - V 2 )V[a\}p + T r V[a 2 }p 
+ Tir?{V[ba\]p + V[tfa\]p) 
+ Tir) 2 {bh[a u a\p] + \pai,a\}b ] b), (18) 
with the notation T> for any operator A: 

V[A]p = ApA* -hA*Ap + pA*A]. (19) 



Here T a = 2-Kp a Vt? a is the rate for electron tunneling to 
the electrode a, while p a denotes the density of states at 
the electrode a. The Liouvillian operator £d describes 
the intrinsic dissipation of the MR induced by its thermal 
bath and can be written in a Lindblad form— as 



£dP 



n{uj m ) + l][2bptf - (tfbp + ptfb)} 



+ ±n(GJ m )[2tfpb-(btfp + pbb% (20) 

where 7 = to m /Q is the intrinsic dissipation rate of the 
MR and n(u) m ) is the average boson number in the ther- 
mal bath. 

Next, we focus on the regime in which the driving 
strength r]Q is low enough so that the time scale related 
to the coupling between the MR and the DQD is slow 
compared to the dynamics of the DQD and the mechan- 
ical oscillation period, just like the Lamb-Dicke regime 
considered in laser cooling of an atom or a trapped ion. 
We also assume that [n(uj m ) + l]j is much smaller than 
the decay rate of the DQD and the oscillation frequency 
of the MR, which is required for appreciable cooling (see 
Sec. IV). In this Lamb-Dicke regime, the DQD can be 
regarded as a structured environment and can be adi- 
abatically eliminated. 48 Since we arc interested in the 
behavior in the limit t — > 00, we can project the system 
on the subspace with zero eigenvalue of £0 (£0 is the 
Liouvillian for the decoupled MR and DQD), according 
to 



Vp = p s d ®Tv d {p} = p s d ®p 1 



1-V, 



(21) 



with p d denoting the stationary (hence the "s" super- 
script) density matrix of the DQD (hence the "d" sub- 
script), and p the density matrix of the MR. Up to sec- 
ond order in 77, the master equation (jTTJ) can be written 



25.48 



dp 
~dt 



£(t)p = [£ (t) + d(t) + C 2 (t)]p, (22) 



(17) where 



(23) 



A>/> = -i[uJ m tfb, p] - i[-Acr z + -~-a x , p] 

+r l V[a\]p + r r V[a 2 }p, 
C lP = C+{t)p + Ci{t)p, 

L\p = -ir]—[(a + -a-)b i , p], 

CiP = iri^-[{<J + -a-)b, p], (24) 

C 2 p = - V 2 T l V[a\}p + v 2 Ti('D[bal}p + -D[tfal}p) 

+rj 2 Yi(b' [ b[ai,a\p\ + \pai,a\]b%) + CbP, 

(25) 

are the Liouvillians to zeroth, first, and second order in 
77, respectively. 



Projecting the master equation in Eq. ([22)) on the V 
subspace, one has 



dt 



Vp = [VC 2 V + Pdi-Co^dVp]. 



(26) 



This result is obtained to second-order perturbation in r\ 
(see Appendix). Since we are interested in the dynamics 
of the MR, we will trace Eq. (|2"6")l over the DQD degrees 
of freedom. The master equation for the reduced density 
matrix of the MR is given by (see Appendix B) 

fi = -i[u) m + 5 m , tfb] + ~{7Kw m ) + l] + A_{uj m )} 
x [2bptf - (6 f 6/x + nb%)] 
+ -{ in {u m ) + A+. (w ro )][2&V& - (W>V + IJM>% 

(27) 

where 5 m is the driving-induced shift of the mechanical 
frequency. 

From Eq. (127|) . one finds that, besides the effects in- 
duced by the coupling to the thermal bath (terms pro- 
portional to 7), the cooling and heating induced by the 
inelastic scattering processes of the MR can occur and 
the corresponding rates A^ are 

A _(oj) = ^°Rc 
v ' 2 

r^o, 



1 



A 2 + ( 7o - w)[( 7 o - ioj) + -fL(w)\ 



x[^-(P(uj)(a y ) s + iR( U )(a x ) s ) 

-(70 - tw)(pj + pI) s + tA(/£ - Pd)s] J + 2D, 

A+(uj) = A_(-w), (28) 

where (• • -) s means the steady state solution of the cor- 
responding quantity and 

_ . . r r + 2r ; - 2iuj 



i^yuj, 


(Ti — iu>)(T r — iu>) ' 




P( W ) = 


7o - iu 
Ti - iuj 


' 2r, (2r ; - ioj) 

—iuj (T r — iuj) 


R(u>) = 


7o -iu 
Ti - ILO 


I ! (2r t -«w)l 

(r r — iuj) 


) 


D = 


^r,(pS 


>- 





(29) 



From Eq. (|2"8)) , one can see that the rates A T depend on 
the tunneling rate Tu r \ . For simplicity, below we consider 
the symmetric couplings of the DQD to the electrodes, 
i.e., Ti = T r = T. 



IV. STEADY-STATE AVERAGE PHONON 
NUMBER IN THE MECHANICAL RESONATOR 

A. Cooling condition 

Below we study the cooling limit regarding the steady- 
state average phonon occupancy in the MR. The equation 
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FIG. 3: (Color online) (a) The rate W/rf as a function of the 
driving detuning A = ujd — luq, for different driving strengths 
Q - For red detuning A < 0, the rate W > and the MR 
can be cooled. In contrast, the heating process dominates for 
blue detuning A > 0. (b) Microwave-induced steady phonon 
occupancy n/ as a function of the normalized driving detuning 
A/r. The positive (negative) value of the steady phonon 
occupancy on the left (right) side indicates that the scattering 
process induced by the microwave produces cooling (heating) 
of the MR. The parameters are r = 1, u m — V, and £7o = 
T (black), Q, = 2 T (red), fio = 5 T (green). 



of motion describing the phonon occupancy distribution 
can be obtained from the master equation (|2"T|) of the 
MR, i.e., 



dp n 
dt 



= {j[n(Lu m ) + 1] + A_}[(n + l)p n +i - np n ] 
+ [yn(uj m ) + A + ][np„_i - (n + l)p n ], (30) 



where p n = (n\fj,\n). At steady state, dp n /dt = 0. Rs 
solution gives the steady-state average phonon occupancy 



(n) 



7n(w m ) + A 4 



(31) 



7 + 1F 
where W = A- — A + is the rate of cooling or heating. 



With the expression (|3Tj) for the steady-state average 
phonon occupancy, we can obtain the cooling condition 
for the single MR. 

To achieve cooling, the condition W > is required. 
Otherwise, cooling the MR is unachievable since the heat- 
ing processes plays a dominant role for W < 0. In 
Fig. EJa), the rate W is plotted as a function of the 
driving detuning for different driving strengths. It can 
be seen that the sign of the rate W exhibits a depen- 
dence on the detuning A. When the driving is red- 
detuned (A < 0), the rate W > and the MR is cooled. 
The steady-state average phonon occupancy is plotted 
in Fig. [3^b). It clearly reveals the cooling (A < 0) and 
heating (A > 0) regions regarding the driving detuning. 
Below we focus on the red-detuned region to discuss the 
cooling of the MR. 

At the beginning, if the MR is in its thermal equilib- 
rium state, the initial average phonon occupancy of the 
MR is given by n(uo m ) — nth- The numerator of Eq. (|3ip 
reveals that two parts contribute to the steady-state av- 
erage phonon occupancy. The first term, jn(uj m ), re- 
sults from the thermal bath and the steady-state average 
phonon occupancy is proportional to the initial thermal 
occupancy. The second term, A + , in Eq. (|3"Tj) originates 
from the scattering process induced by the driving mi- 
crowave. In order to achieve an appreciable cooling, i.e., 
(n) <C n t h, we need W 3> 7. In this regime, the steady- 
state average phonon occupancy is approximately given 
by 

W 



(n) « nf = 



(32) 



In Fig. [5{a), we plot n/ as a function of both the driv- 
ing detuning and the oscillation frequency of the MR. 
The region for cooling the MR is enclosed by the contour 
line rij = 1 . This region covers a wide area in the uj m - A 
plane, implying that the cooling of the MR is experimen- 
tally accessible. Moreover, the region with steady-state 
average phonon occupancy much smaller than unity (e.g., 
?i/ = 0.01 <C 1) is explicitly shown in Fig. rjjb). This re- 
gion corresponds to the cooling of the MR to the ground 
state and it can be achieved by both changing the decay 
rate of the DQD and detuning the transition frequency 
of the DQD from the microwave frequency. 



B. Resolved-sideband cooling 

Next, we analytically study the phonon occupancy 
in the resolved-sideband cooling region, i.e., Lu m ^ T. 
In this regime, the MR's motional sidebands are well 
resolved since the natural lincwidth T of the absorp- 
tion sidebands for different mechanical modes are weakly 
overlapping. This enables highly targeted cooling with 
only one mechanical mode. 

The previous semiclassical analyse o 18 ' 19 i 20 ' 37 i 39 show 
that the resonant system (here, the DQD) cannot re- 
spond instantaneously to the mechanical motion. Hence, 



100 
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FIG. 4: (Color online) (a) Microwave-induced steady-state 
average phonon occupancy nf, shown by colors, as a func- 
tion of both the normalized driving detuning A/r and the 
normalized oscillation frequency uo m /T. (b) Three contour 
curves of the microwave-induced steady-state average phonon 
occupancy for n/ = 1, 0.1 and 0.01. The parameters are 
r = 1, and Q = 2 T. 



this finite response time induces a phase lag, which pro- 
duces a force opposing the mechanical motion, leading to 
a reduction of the mechanical motion. Here, in quantum 
theory, we will see how the energy is exchanged between 
the DQD and the MR during the cooling process (see also 
Ref . |28() . We find below that when the DQD is tuned, via 
varying the gate voltage, to satisfy the lowest sideband 
condition, i.e., ui m = ujQ — uJd, the MR can be cooled suffi- 
ciently. In this case, the anti-Stokes process is resonantly 
enhanced, as discussed in the introduction. 

In the resolved-sideband cooling regime, the rates Azp 
are approximately given by 



iM Wm ) » v 2 r(p° d ) s + 



ffng 



(Aa-c^ + figr + ^iy 



- (a 2 + «*, + n§) {p 2 d + P \) s T Aw ro r( p i - P 2 d ) s 



(33) 
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A = — u) m , the anti-Stokes process (|l)|n) — > |2}|n — 1}) is 
resonantly enhanced. Therefore, an appreciable cooling 
is achieved. The result obtained here is also consistent 
with previous theoretical predictions for the resolved- 
sideband cooling limit ) 25 ! 44 which has been verified by 
experiments! 29 ' 50 However, an important advantage in 
the present set-up is that both the decay rate and the 
energy splitting of the DQD can be tuned by varying 
the gate voltage to reach the resolved-sideband cooling 
regime. In this regime, the steady-state average phonon 
occupancy is much smaller than unity. Thus, the MR 
can be cooled to its ground state. 

Moreover, in Fig. \5\ we plot the steady-state average 
phonon occupancy in the appreciable cooling regime us- 
ing the exact [Eq. (|28|) ] and the approximate expressions 
[Eq. ([55)) ] for the rates A T , respectively. In the resolved- 
sideband regime, the results agree well with each other 
in these two cases [Fig. H[a)]. In Fig. G3h), when the 
coupled MR-DQD system deviates from the resolved- 
sideband regime, however, the exact and the approximate 
expressions can differ significantly from each other. This 
regime implies the breakdown of the approximation used 
to obtain Eq. (|3"3")) . This indicates that one needs to use 
the exact expression (j28[) to describe the MR cooling in 
the non-resolved- sideband regime. Indeed, as shown in 
Fig. [4] though the condition u> m ^ T is not fulfilled in 
the non-resolved-sideband regime, the ground-state cool- 
ing of the MR is still achievable. 



FIG. 5: (Color online) Comparison of the steady-state aver- 
age phonon occupancy n/ between the exact (|28[) and approx- 
imate (|33[) expressions as a function of the normalized driving 
detuning A/cu m . The parameters are r = 1, fio — 2 T, as well 
as (a) LOm/T — 100, and (b) w m /r = 1. 



and the steady-state average phonon occupancy becomes 
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(34) 



As shown in Eq. ([34)l . the steady-state average phonon 
occupancy n/ only depends on the normalized detuning 
A/T for a fixed MR oscillation frequency. With an opti- 
mal detuning value A = — u> m , the minimum limit of n/ 
is found to be 



7 / r 
n min = min{n/} = - 1 

8 V u m 



(35) 



which is much smaller than unity in the resolved- 
sideband cooling region, i.e., ui m 3> Y. This result also 
shows that when the DQD is tuned to the lowest sideband 



C. Estimates 

Finally, let us estimate the steady-state average 
phonon occupancy of the MR in the resolved-sideband 
regime using typical experimental parameters! 29 ' 41 ' 42 
Here we use uj m = 2tt x 100 MHz, r = 2tt x 10 MHz, 
O = 2tt x 20 MHz, and n = 0.2. When the DQD is 
tuned to the lowest sideband A = —u> m , the microwave- 
induced steady-state average phonon occupancy is given 
by 



71/ 



and the cooling rate is 



7/ r 



W = A_ - A_, 



0.00875, 



1.4 MHz. 



(36) 



(37) 



Considering a MR with a quality factor Q = 10 5 , the 
intrinsic dissipation rate of the MR is 7 = 1 kHz. Hence, 
an appreciable cooling effect, i.e., W ^> 7, can be pro- 
duced. For a MR precooled by a dilution refrigerator 
to a temperature T of, e.g., 100 mK, we have nth ~ 21, 
and thus the steady-state average phonon occupancy part 
that comes from thermal fluctuations is 



7 "-th 
W 



0.015. 



(38) 



It follows from Eq. (|31j) that the steady-state average 
phonon occupancy of the MR is given by 



(n) « 0.024 < 1, 



(39) 



which corresponds to an effective temperature T e g rs 
1.3 mK. This means that the cooling of the MR to the 
ground state is achievable using the proposed set-up with 
a MR coupled to a DQD. 

As derived in Ref. |28l , starting from an initial tempera- 
ture To, the final temperature of the cooled MR is bound 
by 



Tf = 



To 



Wo 



(40) 



i.e., this expression provides the lower limit of the tem- 
perature that can be achieved via the sideband cooling. 
Experimentally, the energy level difference Ae between 
the ground state and the first excited state of a single 
quantum dot can be ~ 250 /neV (see, e.g., Ref.l5ll). which 
corresponds to a frequency v ~ 60 GHz. This level dif- 
ference can be even larger by decreasing the size of the 
dot. In our set-up, Ae should be larger than the en- 
ergy level difference Tiloq between the ground state of the 
left dot and that of the right dot, so as to prevent the 
electron in the right dot from tunneling to the first ex- 
cited state of the left dot. Here, for example, we can 
choose ojq = 2ir x 40 GHz. Other parameters used for 
calculating (n) in Eq. (f3"9"| are oj m — 2ir x 100 MHz, and 
T = 100 mK. From Eq. gDJ), we have T f = 0.25 mK. Ob- 
viously, the temperature limit is lower than the achieved 
temperature T c g fw 1.3 mK in Eq. (|39p . This implies that 
the MR can be further cooled using our sideband cooling 
proposal. 



V. CONCLUSION 

We have studied the cooling of a MR by electrostati- 
cally coupling it to a semiconductor DQD. Here the DQD 
works as a two-level system and the decay rate of the 
DQD corresponds to the rate of the electron in the DQD 
tunneling to the electrode. This tunneling rate and the 
energy level splitting of the DQD can be tuned by vary- 
ing the gate voltage. We show that when the two-level 
system is driven by a microwave field in red-detuning, 
the MR can be cooled, in analogy to the laser sideband- 
cooling of atoms or trapped ions in atomic physics. Also, 
we obtain analytical results for the resolved-sideband 
cooling of the MR. Moreover, our results show that the 
ground-state cooling of the MR can be achieved both by 
detuning the transition frequency of the DQD from the 
microwave frequency and by changing the decay rate of 
the DQD. Importantly, this frequency detuning and the 
decay rate of the DQD are tunable by varying the gate 
voltages of the DQD. Thus, the coupled MR-DQD sys- 
tem provides an experimentally implementable set-up for 
ground-state cooling of MRs. 
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APPENDIX A: MASTER EQUATION FOR THE 

DENSITY MATRIX PROJECTED ON A 

SUBSPACE OF THE LIOUVILLIAN Co 

In this appendix, we derive Eq. (|26p by projecting the 
master equation (|17[) on the subspace, with zero eigen- 
value Ao = 0, of the Liouvillian Co for the decoupled 
resonator and DQD system. The projection is defined as 

-pp = fi®Tr d {p} = (f d ®H, Q = l-V, (Al) 

where V is the projection operator. The definition of V 
implies that CqP = V Cq = 0, which leads to 

VC o V = QC o r = VCoQ = 1 QC Q = C . (A2) 

With these relations, the projection of the master equa- 
tion (fTT)) gives 

Vp = PC 2 Vp + P(C 1 +C 2 )Qp, (A3) 

Qp = [C + Q(C 1 + C 2 )}Qp+Q(C 1 + C 2 )Pp. 

(A4) 

Here we have used the relation VC\V = 0, due to the 
fact that tracing over the interaction between the DQD 
and the MR equals zero. Next, we define 



v(t)=Vp(t), w{t) = Qp(t). 
Applying the Laplace transform 

b(s) = J b(t)e- st dt 
" 



(A5) 



(A6) 



on Eq. (|A4[) . one has 

sv(s) - v(0) = VC 2 v{s) + V{d + C 2 )w(s), 
sw(s) - w(0) = [C Q + Q{d + C 2 )]w{s) 

+ Q(£i + C 2 )v{s). (A7) 

Then we introduce a small parameter £ to characterize 
the order of the Liouvillians^ i.e., 

Co (t) -> Co (t) , & - £Ci (t) , C 2 - C 2 C 2 (t) . (A8) 



Substituting Eq. (|A8|) into Eq. (|A7|1 and including terms 
up to second-order, one has 

sv(s) - [v(0) + V{Cd + Q 2 C 2 ){s - £o)- l w(0)} 
= C 2 PC 2 v(s) + fTdis - Co^QCivia). (A9) 

Finally, neglecting the correction due to the initial con- 
dition and performing the inverse Laplace transform, one 
obtains 

Vp = VC 2 T P + VC l Q{-Cu)- 1 QC 1 Vp, (A10) 

which is just Eq. ([26]) in Sec. III. 



and the rates A T (uj m ) 



APPENDIX B: MASTER EQUATION FOR THE 

REDUCED DENSITY MATRIX OF THE 

MECHANICAL RESONATOR 

Below we derive the master equation for the reduced 
density matrix p of the MR. Following the procedures in 
Ref. |48|, we trace over the DQD degrees of freedom in 
Eq. (jAlOj) . the first term on the r.h.s. of Eq. (|A10|) gives 



Tr, 



t {VC 2 V) = -i^°((r x ) a [6t 6, M ] 



+ ^{l[n{u Jm ) + l]+f 1 2 T l (p d ) s } 

x[2bptf -{tfbp + ptfb)] 

+i[ 7 n(u/ m ) + rf > r l (p d ) a ][2tfnb - (6&V + pbtf)]. 



(Bl) 



The second term on the r.h.s. of Eq. (IA10|) gives 



Tr d {VC 1 (-£o)- 1 £iV P }= j dtTtdtP&e-^CiPp} 

o 
= -i lm[S(oj m ) + S(-uj m )][tfb,fj,] 
+ Re[S(w TO )][2tyj&t - (tfbp + ptfb)} 
+ Re[S(-uj m )}[2b^ib - (6&V + nbtf)], (B2) 

where 

+00 

SM^ 2 ! f dte^(* y (t)a y (0)). (B3) 



Substituting Eqs. ((Bl]) . (|B2]) . and |fMJ into ifATO]) . we 
obtain the master equation for the MR 

(i = -i[u m + 5 m , tfb] + - {y[n{u m ) + 1] + A_(w m )} 



x [261*6+ -(tfbu + fitfb)] 

+-[jn(uj m ) + A+{uj m )][2tf pb - (6&V + /*66 f )]. 



This is just Eq. ([27]) in Sec. III. 
In Eq. (|B4[) . we have introduced 



(B4) 



2 fio 



$ m = r) —{<j x )s + Im[S'(w m ) + S{-u m )] (B5) 



with 



A T (w m ) = 2Re[5(±D 



£>=2^(p2>a 



D], 



(B6) 



(B7) 



where (p") s is the probability of an empty DQD at the 
steady state. 

To determine the rates A T , we need to calculate the 
correlation function S(u>) of Eq. (|B3j) using the equation 
of motion for the DQD: 

p d = -i[Aa z +—a x , / 9 d ]+r2£>[a|] j o rf + r,.P[a 2 ]pd- (B8) 



From Eq. (|B8|) . one can obtain the following equations of 
motion: 

'l t-i "n 1 "n 2 / \ 

Pd = T i- T iPd- T Wd-^{vy), 
Pd = - T rP d + Y\ a v)> 

(a x ) = A(a y ) -7o((Tx), 

(a y ) = -A(a x )-n o (p d -p d )- l0 (a y ). (B9) 

where 70 = T r /2. With the normalization condition p d + 
p d + Pd = 1) the steady-state solution of Eq. (|B8|) is given 
by 



2 F 



M ' 

2r/r r iio 

M : 



(P d )s 

(0- X )s 



4r z A 2 + r,ng + r^r 2 



M 



41^0 
M : 



(<Ty)s 

where 

m = 4r,A 2 + n 2 Q (2Ti + r r ) + r ; r 2 . 

From Eq. (|B9[) . one can obtain 



(BIO) 



(Bll) 



(a y (s))= A 2 + ( s + 7 ) 2 + ^( s + 7 ) 



nL s 2s + 2r ( + r r 



xln t 



2 v ~ ' "(s + Ti)(s + T r ) 
r K s + 7) , n (s+7) 1 



s(« + r,) ' °°(s+ro 



:pS(o) 



-n 



(s + 7)(s + 2ri) 



p 2 M 



(« + r,)(s + r r ) 

+(s + 1 )(a v )(0)-A(a x )(0) 



(B12) 



which is the Laplace transform of (a y (t)). Then, using 
the quantum regression theorem,— and performing the 
inverse Laplace transform, one can calculate the correla- 
tion function in Eq. (|B3|1 . Subsequently, with Eq. (|B6|) . 
the rates A-^(u> m ) are obtained as in Eq. 
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